Ground state cooling is not possible given initial system-thermal bath factorization 
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In the quantum regime, ground-state cooling of a small 
object that is embedded in a thermal environment, such 
as neutral atoms [l| , ion traps @ , mechanical resonators 
[1, H|, nuclear spins (polarization) Q, is an intriguing 
challenge and one of the most desirable of quantum tech- 
nologies 0. Mathematically, a ground-state cooling (or 
polarization) process can be formulated as a transforma- 
tion from the initial state of the system+bath to a final 
state, where the small object, the "system", reaches its 
ground state. There have been variety of ground state 
cooling schemes, for example, sideband cooling which 
have been carried out experimentally [7|, |8j. Here we 
prove that a fundamental constraint on the cooling dy- 
namic implies that it is impossible to cool, via a unitary 
system-bath quantum evolution, a system that is embed- 
ded in a thermal environment down to its ground state, 
if the initial state is a factorized product of system and 
bath states. The latter is a crucial but artificial assump- 
tion often included in many descriptions of system-bath 
dynamics @. The analogous conclusion holds for "cool- 
ing" to any pure state of the system. 

To prove this fundamental statement we consider a 
generic arrangement with a small entity, comprising a 
few degrees of freedom, referred to as the "system" pos- 
sibly subjected to time dependent fields, interacting with 
a bath that is in a thermal equilibrium state. The total 
Hamiltonian is given by 



H = H S + H E 



H 



SB, 



(1) 



where H$ is the Hamiltonian of the system, Hb is that 
of the thermal bath, and H$b denotes the system-bath 
interaction Hamiltonian. The details of these terms, 
whether controllable or uncontrollable (time-dependent 
or not), do not alter our results. Time evolution in 
quantum mechanics is dictated by a propagator U(tf,t ), 
which transfers the full initial system+bath density ma- 
trix p(to) to the final density matrix p(tf), 



p(t f ) = U(tf,t )p(t )U^(tf,t ) 



(2) 



Since H is Hermitian, the unitary condition 
U{to,tf)W(t f ,t ) = W(t f ,t )U(t f ,t ) = I is sat- 
isfied, I is the unity operator, and the trace of the 
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density matrix p is preserved. Without loss of generality, 
the initial state of the total system is assumed here to be 
diagonal, p(to) = diag(P , Pi, ■■■)■ If it is non-diagonal, 
we can diagonalize it by a unitary operator W such that 
diag(P , Pi, ...) — Wp W^. The set {P} corresponds to 
the eigenvalues of p(to), and we order the eigenvalues ac- 
cording to their magnitude, 1 > Pq ^ Pi ^ .... Similarly, 
without loss of generality, we can also assume a diag- 
onal form for the final state, p(tf) — diag(Qo 7 Qi, ■■■), 
again ordered as 1 ^ Qq ^ Qi ^ .... The set {Q} 
corresponds to the eigenspectra of p{tf). If p(tf) is not 
diagonal, it can be diagonalized with a unitary matrix, 
diag(Q , Qi, ...) = V^pft^V. Overall, we can redefine 
the time evolution operator U as VUW , a unitary 
operator, to ensure that both initial state and final state 
are diagonal matrices. The density matrices p(tf) and 
p(to), Hermitian operators, are connected by a unitary 
(rotation) operation, thus they must have identical 
eigenspectra, i.e., {Q} = {P}. Since the elements are 
ordered, we can relate them one by one, 



P k = Qk, k = 0,1, 2. 



(3) 



We now define do and df as the number of non-zero eigen- 
values in the set {P} and {Q}, respectively. Under a uni- 
tary evolution, do must be equal to df, a prerequisite for 
relation to hold (or do/df = 1 as do goes to infinity). 

Based on these simple considerations, we argue next 
that under system-bath unitary operations, acting on 
system or bath or both, one cannot cool a mixed sys- 
tem state down to its ground state if the total density 
matrix is initially system-bath factorizable and the bath 
is thermal. That is, the process 



p{to) 



b 



B = p{t f ) 



(4) 



cannot be carried out with a unitary matrix U (even if 
it operates on both the system and bath) . The left hand 
side in Eq. describes the initial system-bath prod- 
uct state. Here, s — diag(sn, si, ■■•) denotes the system 
density matrix at to, which is anything but a pure state, 
and b = diag(&o> b±, ...) denotes the bath state at that 
time, a thermal state at nonzero temperature. The right 
hand side of Eq. ([4| includes the target final state where 
the system has been cooled down to its ground state |0) 
and the bath is a mixed state B = diag(i?o, B\, ...) which 
is a diagonal matrix that does not necessarily describe 
a thermal state. We now provide an argument, which 
shows that one cannot evolve between these initial and 
final states via unitary dynamics. 
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Define N$ and Nb as the Hilbert space dimension of 
the system and bath, respectively. If the bath is initially 
thermal and s is not a pure state, the inequality do > Nb 
holds. In particular, if the system is initially thermalized 
we reach the upper bound do — N$Nb- On the other 
hand, the target state p(tf) has only df < Nb nonzero 
eigenvalues, reaching the bound df = Nb if the bath be- 
comes a thermal state at time tf. Since df < do, equation 
Q, written here in the form 

s m bj = B k (5) 

cannot be satisfied. Here the index m counts the system 
eigenvalues, j and k follow the bath eigenvalues. Hence, 
system-bath unitary operations cannot cool a system con- 
nected to a thermal bath if the system-bath initial state is 
factorizable and the system is initially in a mixed state. 
The analogous proof holds for any final pure state of the 
system. 

This proven conclusion stands in sharp contrast to 
a multitude of studies, based on master equation ap- 
proaches, that demonstrate ground state coo ling from an 
initial product system-bath state (e.g., Ref. [Hi). While 
previous studies may have pointed out the unattainabil- 
ity of the absolute zero in such situations 

EMI, 

here 

we isolate the centrality of the factorization assumption, 
and emphasize its strong implications regarding both the 
underlying physics and the reliability of master-equation 
type computational frameworks that often assume fac- 
torization. 

It is of interest to examine a few related situations. 
First, if the system is prepared in a pure state, we find 
that do = Nb, and ground state cooling can potentially 
be performed if Eq. ([5]) is satisfied. Second, one can 
achieve ground state cooling by preparing the bath in a 
non-thermal state. In this case we consider an initial bath 
state b with N' B nonzero eigenvalues, N' B < Nb- This re- 
sults in d > N' B while df < N B - These values could be 
made identical if the states b and B are very different. As 
the simplest example, consider both the system and the 
"bath" as single qubits, where initially the "bath" pop- 
ulates its ground state, bo — 1. Using Eq. ([5]), matching 
eigenvalues, we require that so = Bo- System ground 
state can therefore be reached here by the swapping op- 
eration. A more involved scenario includes a two-qubit 
bath and a single-qubit system where we initially set the 
system in a mixed state while we prepare the bath in 
a non-thermal state with precisely two zero eigenvalues, 
b 2 = and 63 = 0. The prerequisite for ground state 
cooling, do = df, could be fulfilled here if at the end of 
the quantum evolution all four bath eigenvalues B^s are 
made nonzero, resulting in do — df = 4. 

Third, we note that system-bath correlated initial 
states [14[ do allow for cooling. We illustrate this pos- 
sibility by modeling the system as a qubit, with ground 
state |0) and excited state |1). We construct the following 
correlated initial state 

p(io) = |0)(0|®&(°) + |l)<l|®&« (6) 



where = diag(6o, b n ) and bW = 

diag(6 n+ i, 6jv b )- As before, the target state is 
P(tf) = |0}(0| ® B. It is easy to confirm that the 
prerequisite for cooling is satisfied, and the number of 
non-zero eigenvalues for the initial and final density 
matrices is identical, do = df = Nb- Furthermore, 
one could pair the eigenvalues one by one, as required 
by Eq. ([3]). For example, we can set the system with 
so = si = h and the bath with = Bk- As a result, the 
reduced density matrix of the bath is the same, initially 
and finally, whereas the reduced density matrix of the 
system at time to is (|0) (0| + |1) (l|)/2, for a given n 
such that b + ... + b n = 1/2. Note also that cooling may 
also be achieved when measurement is involved in the 
cooling process [l5j . 

Finally, consider approximate ground state cooling, de- 
fined as the evolution from the initial state p(to) — s ® b 
to the final-factorizable state p(tf) — S (g> B, where 
S = diag(So, Si, ...), the diagonal state of the system 
at the final time, describes a system "colder" than the 
initial one, for example in the sense that S has fewer 
nonzero elements than s. This situation is typically as- 
sumed in the framework of Markovian master equations 
[l6| . Since the underlying quantum dynamics is unitary, 
we should be able to match the eigenvalues of the ini- 
tial state and the final state. In particular, the first two 
eigenvalues should fulfill so&o = SqBq and sob\ = S0B1. 
The second relation holds as we assume that the system 
energy gap, between its ground state and first excited 
state, is larger than the corresponding gap in the bath, 

sobi > s\bo- These relations yield B\/ Bo = b\/bo, trans- 
lated to e -(Ei-E )/K BTi = e -( El -E )/K B T f ig Bdtz _ 

mann constant), if we further demand that the bath in- 
ternal spectra is identical at to and tf, and that the bath 
acquires a thermal equilibrium state at the final time. 
Here Ti and Tf denote the temperature at the different 
times. The last relation implies that the final-time tem- 
perature is equal to the initial-time temperature, i.e., the 
bath has not been changed through the cooling process, 
{&} = {B}. As a result, to satisfy Eq. ©, we must con- 
clude that the system retains all its values, S m — s m . In 
the scenario described here, quantum evolution cannot 
modify the system population. Thus, even an approxi- 
mate cooling is impossible, as long as the system ground 
state is nondegenerate. 

In summary, ground state cooling within system-bath 
unitary operations is not possible given initial system- 
thermal bath factorization. The linearity of unitary op- 
erations has, in the past, resulted in a no-go theorem, 
the no-cloning theorem [17| . one of the building blocks 
in modern quantum information theory. Our no-go prin- 
ciple is similarly based on unitary evolution, and it lays 
down the foundation for any theory that aims at describ- 
ing ground state cooling and pure state preparation. For 
example, many recognized master equation techniques, 
as well as Kraus operator based methods @, are predi- 
cated on the initial factorization of the system and bath. 
Adopting these approaches to address issues of cooling 
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should be done with extreme caution, considering the 
fundamental constraint exposed in this work. 
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